MPP-2011-21 



Supernova bound on keV-mass sterile neutrinos reexamined 
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Active-sterile neutrino mixing is strongly constrained for rris > 100 keV to avoid excessive energy 
losses from supernova cores. For smaller m^, matter effects suppress the effective mixing angle 
except for a resonant range of energies where it is enhanced. We study the case of Ut-Vs mixing 
where a i/t-i't asymmetry builds up due to the strong excess of Vs over i/s emission or vice versa, 
reducing the overall emission rate. In the warm dark matter range < 10 keV the mixing angle 
is essentially unconstrained. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Sterile neutrinos i^s can be produced in the early uni- 
verse or in supernova (SN) cores if they mix with one of 
the active flavors P, . Even if the mixing angle 6 is very 
small, repeated collisions of the active component allow 
for an efficient Vs production. The SN 1987A neutrino 
signal duration implies sin^ 29 < 10~^ to avoid exces- 
sive energy losses [3 10]. In detail, this limit depends on 
whether the dominant mixing is with v,, or one of the 
other active flavors. This result assumes rus ^ 100 keV 
for matter effects on neutrino propagation to be negligi- 
ble compared with the vacuum mass. 

For smaller masses, the matter effect typically sup- 
presses the effective mixing angle and thus diminishes 
the limit on sin^ 29. For rUs ^ 1 keV even maximal mix- 
ing is allowed. However, no detailed treatment of the SN 
bound exists for 1 keV ^ to^ ^ 100 keV, where mat- 
ter and resonance effects are important [j, @. On the 
other hand, this is precisely the mass range where ster- 
ile neutrinos coul d play a n interesting warm dark matter 
role in cosmology [lll - [l3| . While sterile neutrinos can be 
produced in the early universe by different mechanisms 
[Til [12} , the production by oscillations and collisions once 
more depends on the active-sterile mixing angle, so one 
naturally wonders about the SN bound on 9. 

In particular, we study the role of feedback of sterile 
neutrino emission on the emission rate itself. When mat- 
ter effects are important, the mixing angle is resonantly 
enhanced in some range of neutrino energies [l3 - [l6j . At 
first the emission of Ds is more efficient than because 
in the neutrino sector the mixing angle is suppressed for 
all energies. As a consequence, the active flavor, being 
trapped in the SN core, builds up a. v excess. The op- 
posite could happen for larger mixing angles, when most 
antineutrinos Ds are trapped rather than freely escape 
and the emission of ly^ is more efficient than that of D^. 
To be specific we use Vt as the active flavor because there 
is no initial Vt-i^t asymmetry in a SN core and TOt- is so 
large that charged t leptons never play any role. Our 
main point is that the depletion of or Dr relative to 
the other always goes in the direction of quenching the 
initial emission rate, implying that in the 1-100 keV-mass 
range the 9 bounds are indeed suppressed. 



Our conclusion differs somewhat from previous discus- 
sions where it was suggested that the buildup of a v-y-v^ 
asymmetry goes in the direction of reducing the matter 
effect and leads to restrictive SN limits on 9 [l3|. While 
we agree that the matter effect can be modified in this 
direction, the positive Vt-Vt asymmetry also implies a 
depletion of source Vr to be converted to relative to 
source to be converted to Vg. It is also possible in 
some region of the parameter space that the matter ef- 
fects are enlarged due to a negative Vt-Vt asymmetry. In 
both cases, the compound effect is a reduction, not an 
enhancement, of the energy loss. 

In Secini we briefly review the matter effects on active- 
sterile neutrino mixing in the SN core. The development 
of a Vr-v-r asymmetry is discussed in Sec. IIIIl where we 
identify a stationary state with equal neutrino and an- 
tincutrino emission rates and estimate the time scale to 
reach it. Sec. |IV] is devoted to the calculation of the 
energy-loss rate caused by sterile neutrinos and the SN 
bound on sterile neutrino masses and mixing angles. Fi- 
nally, we summarize our conclusions in Sec. [V] 



II. MATTER EFFECTS 

The dispersion relation of neutrinos will be modified 
in matter due to the coherent forward scattering of neu- 
trinos off background particles [l^. This matter effect 
can be described by an effective potential Vv^ for each 
kind of active neutrino = Ve, ^'/^ and Ut- The effective 
potentials for antineutrinos have the opposite signs, i.e., 
Vj}^ = — Vj/^. In the case of Ut-Vs oscillation in matter, 
the effective Hamiltonian is 



HcS — 



Yy^ - UJC20 



(1) 



where S29 = sin 26', C20 = cos 29 with 9 being the vacuum 
mixing angle and w = Am^/2i? the oscillation frequency 
in vacuum. As far as the keV-mass sterile neutrinos are 



concerned, we have Ar 



with rris being the sterile 



neutrino mass. In contrast with the flavor conversions 
of ordinary neutrinos in the SN, there are no collective 
effects in active-sterile neutrino oscillations. 
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In SN cores, the main ingredients of matter are pro- 
tons p, neutrons n, electrons e, some muons /i as well as 
active neutrinos i>a and antineutrinos 9a- We ignore the 
possibility of a meson condensate or hyperons. In such a 
medium, the effective potential for tau neutrinos is 



V2GfNb 



1, 



(2) 



where Gf is the Fermi constant, A^b the baryon number 
density, and Yx = {Nx - Nx)/Nb with Nx and Nx 
being the number densities of particle X and its antipar- 
ticle X. Because of charge neutrality Yp ^ Yg + F^, 
we have F„ = 1 — — 1^. While ly^ and of 
all flavors can be produced in pairs, electron neutri- 
nos i/e can also be generated in beta processes such as 
e~ + p ^ Ve + n, and similar for muon neutrinos. The 
muon mass = 106 MeV is comparable to the average 
thermal energy (E) = 3T with T = 30 MeV, thus a small 
population of muons is unavoidable. Beta equilibrium 
j/^ + n o p^ + p implies the relation among chemical 
potentials p^ — fi^^^ = pLn — pLp = jl, where /t typically 
lies in the range 50-100 MeV, depending sensitively on 
the equation of state. Noting that the initial /i lepton 
number is vanishing and taking /i = 50 MeV, one finds 
Pp « 18 MeV and p^^ « -32 MeV On the other 

hand, the v,, chemical potential will be much larger due 
to the electron lepton number trapped during infall. 

For simplicity, we consider a SN core just after the 
bounce and assume the temperature T = 30 MeV and 
matter density p = 3.0 x 10^^ g cm~'^ to be constant. 
We furthermore take a typical value of electron lepton 
number fraction Yl ^ Y^, + YJy^ ~ 0.37, which leads to 
Ye = 0.3 and Y,y^ — 0.07 due to beta equilibrium with 
/t = 50 MeV. It is straightforward to verify that Yf^/Y^ = 
0.01 and Y,^^/Y,^^ = —0.05, so henceforth we simply set 
Yp = Yi, ^ =0. Tau neutrinos initially follow the Fermi- 
Dirac distribution without chemical potential, but later 
an asymmetry develops due to Vs emission. With these 
simplifications, the effective potential in Eq. 1^ is 



= -^^B (1 - Ye - 2Y,^ - 4y,J 



(3) 



It is now evident that V^^ is negative for Y^,^ — 0, imply- 
ing that the Mikheyev-Smirnov-Wolfenstein resonance 
occurs in the antineutrino sector ^lA . .15j . 

Given the effective Hamiltonian in Eq. ([1]), one imme- 
diately obtains the effective mixing angle 



sin"' 29,, 



sin^ 29 



sin^ 29 + {cos29 ± E / Erf ' 



(4) 



where the upper sign refers to v and the lower to v. The 
resonant energy E^ = Am? /2\V,y^ \ is 



E, = 3.25 MeV 



10 keV 



p^l\Y^-Y,.^ 



(5) 



where pi4 is the matter density p in units of 10^"* g cm ^ 
and Yo = {l-Ye- 2y^J/4. As indicated by Eq. (g]), the 



mixing angle 9^ for the whole energy range is always sup- 
pressed by matter effects, while 0p can be resonantly en- 
hanced for E ^ Er cos 29. Note that the "vacuum limit" 
with 9y K 9ij K 9 is reached for large sterile neutrino 
masses ^ 10 keV, while the "medium limit" with 
reduced mixing angle w 0p w {Er/E)9 is obtained for 
small masses TOs <C 10 keV. For intermediate masses, we 
have a resonance in the antineutrino sector. 

Sterile neutrinos are produced in the SN core via os- 
cillations and collisions of tau neutrinos. If the effective 
mixing angles 9i, and 0p are small enough, and 
can escape from the core immediately after production. 
Since the mixing angle of antineutrinos is always larger 
than that of neutrinos, the emission rate of antineutrinos 
exceeds that of neutrinos. Consequently, a Vt-Vt asym- 
metry arises from these different emission rates. As we 
can observe from Eq. ([5|), it might turn out that a rela- 
tively large Vt-v^ asymmetry is achieved, i.e., Yy^ — >■ Iq, 
so as to drive the resonant energy E^ to infinity, leading 
to the "vacuum limit" even for small masses. Put another 
way, the Vr-^r asymmetry seems to develop in the direc- 
tion of reducing the matter effects, resulting in restrictive 
SN limits on the vacuum mixing angle 9 for both large 
and small sterile neutrino masses. It has been argued [lOf 
that there exists a stationary state with Y^^ = Iq, which 
can be achieved rapidly and thus validates the "vacuum 
limit" and restrictive bounds on 9. 

However, the state with Y^^ — Yq cannot be station- 
ary, because the population of Vt at this moment is larger 
than that of Dr, implying that more neutrinos than an- 
tineutrinos are ready to be emitted and thus the con- 
dition y,y^ = Iq breaks down. Furthermore, it is even 
possible that the emission rate of Vg exceeds that of Pg, 
since the effective mixing angle 0p can be so large that 
most of Vs are trapped in the core. Therefore one may 
obtain a negative asymmetry Y^^ < 0, driving the system 
towards the "medium limit." In the following sections, 
we shall examine how the Vt-Dt asymmetry actually de- 
velops, and explore its implications on the anomalous 
energy-loss rate of the SN core and thus the SN bounds 
on sterile neutrinos. 



III. STATIONARY STATE 

A. Weak-damping limit 

The matter density of SN cores is so high that both 
neutrino oscillations and frequent collisions with back- 
ground particles arc important. An elegant method to 
treat neutrino fiavor conversions in this case is to im- 
plement the matrix of occupation numbers {pp)ij = 
(6j(p)6i(p)), where fei(p) denotes the annihilation oper- 
ator for a neutrino of flavor i and momentum p, and to 
derive the non-Abelian Boltzmann equations of pp (isl - 
[20j . The diagonal elements (pp)ii are the usual occu- 
pation numbers /p, while the off-diagonal ones encode 
the phase information. The analogous definition for an- 
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tineutrinos is {pp)ij = {dl{p)dj{p)) with di{p) being the 
annihilation operator for an antineutrino of flavor i and 
momentum p. In general, the equations of motion for pp 
and Pp are complicated by the nonlinear nature of the 
collision integrals. 

For keV-mass sterile neutrinos, the problem can be 
much simplified by taking the weak-damping limit, which 
is usually valid in SN cores. To be more explicit, we es- 
timate the neutrino oscillation length in matter 



Aosc < 0.7 



E 



^30 MeV J \sm20 J \ m 
and the mean free path of tau neutrinos 



A..p = ^«l.lxl03cm(^)%-, (7) 

where a^M ^ GpE'^/n is the cross section of neutrino- 
nucleon scattering via the neutral-current interaction. 
The weak-damping limit is Aqsc ^ A^fp, meaning that 
the active-sterile neutrino oscillations take place many 
times before a subsequent collision of active neutrinos 
with the nucleons. It deserves mention that the weak- 
damping limit is violated for smaller masses <^ 1 keV, 
however, the effective neutrino mixing angle in this case is 
highly suppressed by matter effects and thus the energy- 
loss rate is negligibly small. We shall always assume the 
weak-damping limit in the mass range of our interest. 



B. Emission rates 

In the weak-damping limit, pp is averaged over many 
cycles of oscillations and then can be parametrized by the 
neutrino occupation numbers and with the neu- 
trino energy E = |p| and similar for antineutrinos [1]. 
Note that we consider a homogenous and isotropic en- 
semble of active neutrinos and antineutrinos, which are 
trapped in the SN core and stay in thermal equilibrium 
with ambient matter, so the momentum direction is ir- 
relevant. Assuming that sterile neutrinos freely escape, 
we obtain the evolution equation of the density 



E^dE 



E'^dE' 
27r2 



W^E'eIe' , (8) 



where a denotes the target particle, and W^,^ the tran- 
sition probability for i>t{E') + a — > h'riE) + a via the 
neutral-current interaction. In a similar way, we can de- 
rive the evolution equation of the Dt number density, in- 
volving the mixing angle 9p, the occupation number 
and the transition probability W^i^- However, only the 
neutrino energies ensuring 9^,0[, ^ 0c are taken under 
the integration, where the critical mixing angle 9c ~ 10~^ 
can be estimated by requiring the mean free path of ster- 
ile neutrinos to equal the core radius R 10 km. On 
the other hand, the transition probabilities W^,^ and 



W%,E of neutrino-nucleon scattering dominate over those 
of neutrino-electron scattering if the sterile neutrino mass 
is lying in the range of 1 keV < ttIs < 100 keV. Hence we 
consider only the neutral-current v-N scattering in the 
leading-order approximation. 



C. Degeneracy parameter 

In order to describe the Vt-v-t asymmetry, we assume 
the occupation number /^(t) — [exp(£'/T — r\(€)) -|- 
for Vr and /^(t) = [exp(_E/r -f r\(t)) + for Vt, where 
T is the temperature of the SN core and r\(t) is the degen- 
eracy parameter. In the absence of active-sterile neutrino 
mixing, the spectra of Vr and Vt just follow the Fermi- 
Dirac distribution with a vanishing chemical potential, 
namely 77 = at i = 0. Substituting the occupation num- 
bers into Eq. (|8]), taking account of the neutral-current 
Vt-N scattering and subtracting the corresponding equa- 
tion for antineutrinos, one arrives at 



At 



An 



- ^Av)W-Hv) , (9) 



where the relevant functions are defined as follows 



dx 



X \ — 3§{x,x-[e ,x-ce^) 

6^+" -f 1 slg + {C20 - X/Xrf 

x^ e{x + x,e-) 
7-^ ; ; — z^dx , (10) 

e^-^ + lslg + {c2e + x/x,)^ 



and ^^(r/) = d[F2{r]) — F2{—ri)]/dr] with i^2('7) being the 
Fermi-Dirac integral of order two. In addition, we have 
introduced x,- = E^/T, = S2(e^±9)/ S2e^, and the box 
function a, b) = Q{x — a) — Q{x — 6), which equals 
one for x € [a, b] and vanishes otherwise. Here Q{x) 
denotes the unit step function, i.e., Q{x) = for x < 
and 8(x) = 1 for a; > 0. Note that Eq. ([10]) has been 
cast into a compact form so as to include both 9 < 9c 
and 9 > 9c cases. 

Taking typical values of the matter density pn = 3.0 
and the core temperature T = 30 MeV, we can rewrite 
Eq. dH) as r]{t) = T^'^J^iri) with ro 1 s (I0"7s2e) and 
J^(n) = Wuiv) - ■^I'irhW'^i'il)- The time evolution 
of the degeneracy parameter ri(t) depends crucially on 
the initial difference between neutrino and antineutrino 
emission rates, i.e., Tq^,M'{0), as well as the evolution of 
J^{ri) with respect to 77. In Fig. 1, we show the initial rate 
77(0) — rQ"^,i^(0) in the (sin^ 26*, ttIs) plane, where the 
blank region with denotes a strong excess of Ds over 
i/s emission while that with "— " represents the opposite 
case. The cyan regions on the left-hand and right-hand 
side indicate 77(0) > and 77(0) < 0, respectively, but 
the magnitude of 1 77(0)! is extremely small for both cases. 
The reason for the former case is just that the effective 
mixing angle is too small, and for the latter case is that 
most of and Ds are trapped. 

If the initial emission rate of Ds is larger than that of i>s, 
i.e., 77(0) > 0, the degeneracy parameter increases from 
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JQ-12 JQ 



10"° 10" 
sin^ (2 6) 



FIG. 1: Density plot of the initial rate 77(0) = Tq^,^{0), 
where the blank region with "+" stands for a strong excess 
of Vs over Va emission while that with "— " is the other way 
round. The narrow strip indicates a rapid and continuous 
transition between these two regions. 



zero to a positive stable value 77*, at which the emission 
rates of Ug and i>s are equal, namely ^^(77*) — ^^{r]*) 
or equivalently Jif{ri*) = 0. In this case, we are finally 
left with a positive Vt-Vt asymmetry when such a sta- 
tionary state is reached. If the initial emission rate of Vg 
exceeds that of Vs, i.e., 77(0) < 0, the degeneracy parame- 
ter decreases from zero to a negative stable value 77* . At 
this moment, we have equal neutrino and antineutrino 
emission rates as well, but a negative Vt-Vt asymmetry. 

The important point here is feedback of the established 
Vt-v^ asymmetry or a finite degeneracy parameter. In 
the case of 77(0) > 0, a positive r\ suppresses the Vr pop- 
ulation and shifts the resonant energy to a larger value. 
The combined result is just to reduce the emission rate. 
Meanwhile, the population of v.^ is accordingly increased 
and the mixing angle becomes less suppressed, enhanc- 
ing the Vs emission rate. Similar arguments apply to the 
case of 77(0) < 0. It is also possible that 77(0) = 0, which 
lies in the narrow strip in Fig. 1, the system remains in 
its initial state with a vanishing Vt-Vt asymmetry. Hence 
we have no feedback effect in this special case. 

In Fig. 2, the degeneracy parameter r\* of the station- 
ary state has been solved from ^vW) = ■^viv): no mat- 
ter whether such a stationary state can be reached. The 
large values of rj* appear in the region where the ini- 
tial emission rate of I's is significantly different from that 
of Dg. This condition can be satisfied for (1) small mix- 
ing angles and intermediate masses, when the mixing an- 
gle 9[> is resonantly enhanced so that =^p(0) » =^i/(0); 
(2) large mixing angles and large masses, when more 
Ds than are trapped in the core such that =^p(0) <C 
^^(0). However, it is obvious that the Vg and Dg emission 
rates in both cases are extremely small. 
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FIG. 2: Contour plot of the asymptotic degeneracy parameter 
77*, which is determined by ^^(77*) = J?^ (?;*), implying equal 
neutrino and antineutrino emission rates. 
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FIG. 3; Contour plot of the estimated timescale r = 
ri* /[tq^ ,j^{0)] to reach the stationary state, where the con- 
tours are labeled by the logarithmic values log(T/l s). 



r-ijr(o) 

= 77V[ro-^^(0)]. 



The timescale for the system to achieve the stationary 
state can be determined by numerically solving Eq. ([9]). 
For a rough estimate, we take r]{t) ^ 77(G) 
and then obtain the timescale r 
In Fig. 3, we show the estimated timescale r in the 
(sin^ 20, rris) plane. In the mass range 1 keV ^ mg < 
10 keV, where sterile neutrinos can be warm dark mat- 
ter, the timescale is larger than the neutrino diffusion 
time Td = 1 s for small mixing angles sin^ 20 < 10"^. 
For larger mixing angles, the stationary state can be 
achieved within a fraction of a second, but rj* in this 
case is quite small as shown in Fig. 2. Hence we ex- 
pect that the feedback effects are negligible. In this con- 
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jQ-8 jQ-6 jQ-4 jQ-2 I 

sin^ (2 6) 

FIG. 4: Supernova bound on sterile neutrino masses rus 
and mixing angles 6, where the purple region is excluded by 
the energy-loss argument while the green one by the energy- 
transfer argument. The excluded region will be extended to 
the dashed (red) line if the build-up of degeneracy parameter 
is ignored, i.e., ri{t) — 0. The dot-dashed (green) line rep- 
resents the sterile neutrinos as dark matter with the correct 



relic abundance Qgh 



0.1. 



nection, the most interesting parameter space should be 
20 keV < TO, < 80 keV and 10"^ < sin^ 20 < lO"'^, 
where both sizable rj* and r < Td are expected. Since 
\J^{t])\ decreases from the initial value \J^{0)\ to zero, as 
the absolute value of the degeneracy parameter |?y(i)| in- 
creases from zero to |77*|, the relaxation time may be un- 
derestimated. However, this rough estimate has already 
shown the main features of the relaxation timescale. 



IV. SUPERNOVA BOUNDS 

The emission of sterile neutrinos and may cause 
rapid energy losses from the SN cores, which can sig- 
nificantly shorten the duration of neutrino signals |2ll |. 
In order to avoid conflict with the observation of SN 
1987A neutrinos, we require the energy-loss rate per 
unit mass to be smaller than 1.0 x 10^^ erg g~^ s~^, 
which can be translated into the volume energy-loss rate 
^ < 3.0 X 10^3 erg cm'^ s'^ for pu = 3.0. Given the 
emission rate of neutrinos in Eq. ([8]) and the counterpart 
for antineutrinos, the energy-loss rate is 

m = mv) + ^.m , (n) 



with 



1 — ^(XjXi-e jXj-e'^) 
e^+n + 1 {c20 - xlx,f 

x^ 6(a;-|-Xre~) 



dx 



where the definitions of relevant parameters are given 
below Eq. ([T0|) . The energy- loss rate S(t) depends on 
time through the degeneracy parameter r](t)^ for which 
the time evolution has been discussed in last section. 

To constrain the sterile neutrino mass and mixing an- 
gle, we evaluate the emission rate 



'd / 

10 



(13) 



averaged over the neutrino diffusion timescale Td = 1 s. 
Beyond the diffusion timescale, one may expect that all 
active neutrinos have already diffused out of the core 
and thus the emission of sterile neutrinos is physically 
meaningless. Our strategy is to follow the time evolu- 
tion of r]{t) for each point in the (sin^2^,TOs) param- 
eter space, and then calculate the averaged energy-loss 
rate in Eq. (IT3l) . We show in Fig. 4 the contour plot 
of the averaged energy-loss rate (t?) in the (sin^26',TOs) 
plane, where the purple region corresponds to {(o) > 
3.0 X lO^'^ erg cm~'^ s~^ and is thus excluded. 

Some comments are in order. First, the initial energy- 
loss rate with 77 = is shown in Fig. 4 for comparison. 
Except for 20 keV ^ to, ^ 100 keV and small mixing 
angles, the averaged energy-loss rate ((f) cannot be dis- 
tinguished from the initial one S'{0). The reason is that 
either the stationary state has not been reached within 
1 s, or the asymptotic value rj* for the stationary state is 
quite small. In the mass range 20 keV < to, < 100 keV, 
the buildup of a i^r-i^T asymmetry is efficient and strongly 
reduces the energy-loss rate. 

Second, keV-mass sterile neutrinos can be produced in 
the early universe and contribute as dark matter to the 
total energy density. In a nonresonant production scheme 
without large primordial lepton asymmetries, the relic 
sterile neutrino abundance can be estimated as [lol [l3j 



0.3 



sin^ 261 



)( 



100 keV 



(14) 



, o / / Noda; , (12) 



The correct dark matter abundance il,^^ = 0.1 is shown 
in Fig. 4, where one can see that the masses around 
TOs — 50 keV have been excluded by the SN bound if 
the degeneracy parameter is assumed to be vanishing. 
In a realistic situation, this region is retrieved because 
the energy- loss rate is reduced as the Vt-Vt asymme- 
try builds up. However, the warm dark matter range 
1 keV ^ TOs ^10 keV is essentially unconstrained. For 
sterile neutrinos of masses below 1 keV, even maximal 
mixing is allowed, because the total energy-loss rate is 
highly suppressed by matter effects. 

Finally, our discussion was based on the energy-loss 
argument, ignoring the sterile neutrinos trapped in the 
SN core. For this reason, the top-right green region in 
Fig. 4, where both the sterile neutrino mass and vacuum 
mixing angle are large, was not excluded. However, the 
mean free path of these trapped sterile neutrinos is much 
larger than that of v^-^ so the energy transfer in the SN 
core will be more efficient. The energy in a SN core is 
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carried by those trapped particles with the largest mean 
free path, a role played here by the sterile states. Once 
more the neutrino burst duration is shortened too much 
if the sterile mfp is larger than a few times the one for 
Vt Q • Therefore, the large mixing angle region is actually 
excluded in the spirit of the energy-transfer argument. 

V. CONCLUSIONS 

Since keV-mass sterile neutrinos are a promising candi- 
date for warm dark matter, we have revisited the super- 
nova bound on the sterile neutrino masses and mixing 
angles by studying the case of Vt-Vs mixing in the SN 
core and requiring no excessive energy losses induced by 
sterile neutrinos. It turns out that the warm dark matter 
range is essentially unconstrained, while sterile neutrinos 
of masses around 50 keV receive the most stringent con- 
straint, i.e., sin^ 2Q < 4.0 x 10^^. For even larger masses 
rus ^ 100 keV, the SN limit on the mixing angles is 
sin^ W < 5.0 X 10~^, which is about one order of magni- 
tude weaker than that for mg ~ 50 keV. It is the matter 
effects that render these constraints quite different. 

We have identified a mass range 20 keV % nis ^ 
100 keV where a sizable v^-Vt asymmetry can be es- 
tablished due to the strong excess of Vg over emis- 
sion or vice versa. The build-up of this asymmetry feeds 
back on the emission rates, leading to a stationary state 
where the neutrino and antineutrino emissions become 



equal. For proper mixing angles, such a stationary state 
can be achieved within the neutrino diffusion timescale 
Td = 1 s. As a consequence, the energy-loss rate will be 
significantly reduced, and thus the bounds are relaxed. 

As for the I'^-i/g- mixing case, our discussions about 
the feedback effects are essentially applicable. However, 
the charged-current interactions of i'^ and should be 
taken into account, and the change of v^-v^ asymme- 
try will be redistributed between muon neutrinos and 
charged muons. The Ue-Vs mixing in SN cores is more in- 
volved because of the large trapped electron number and 
high Ve degeneracy. Besides energy loss, deleptonization 
by sterile neutrino emission is an effect to be taken into 
account. This case requires a dedicated investigation. 

Wc have performed a "single zone" analysis by assum- 
ing a homogenous and isotropic SN core with constant 
matter density and temperature. This treatment should 
capture the dominant feedback effect. However, the local 
variation of these quantities may modify the final results, 
for example smearing out the resonance in the energy-loss 
rate. Such a refinement is also left for future works. 
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